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THE INTEGRAL AS THE LIMIT OF A SUM, AND A THEOREM 

OF DUHAMEL'S 

By William F. Osgood 

1. Introduction. A satisfactory treatment of the integral as tlie limit 
of a sum is to be found in any sbindard text-book on the calculus. If/'(x) 
is a real function of the real variable x, continuous throughout the interval 
a ^ X ^ b, then the sum 

f(x[)Axi +/(xOAa;2 + • • • +f(x'„)Ax„, 

where AXf = Xi — Xj_i, Xq= a, x„ = b, and cc,_i ^ xl ^ cc,-, converges towai'd 
a limit when n becomes infinite, the longest Axi converging at the same time 
toward zero. Furthermore this limit, the definite integral oif{x) taken be- 
tween the limits a and b : 

/ f(x)dx, 
may be computed in terms of the indefinite integral 

F{x) =ff(x)dx 
by means of the theorem that 

''/{x)dx= F{b) - F{a). 



Ja 



This and the corresponding theorem for multiple integrals, I say, are 
treated satisfactorily in stsxndard treatises on the calculus.* The treatment of 
the application of these theorems to problems in physics and mechanics, both 
in courses and in text-books on the calculus, is however far from satisfactory, 
the attitude of the mathematician often being that these applications are without 
value for mathematics, while for the physicist any reasoning is good enough 
which in the long run leads to the right formula. To throw rigor to the 
winds as soon as the hypnotic influence of the cnvh'onment of a course in 



•For an excellent treatment of double and multiple integrals cf. Goursat, Cours d'anali'se 
mathematique, Paris, 1902, ch. G. 

(161) 



162 OSGOOD [July 

pure mathematics ceases is to regard rigor as a frill or as a luxury, as some- 
thing having at most aesthetic reasons for existence, — not as a habit of thought, 
of practical value in research work. On the other hand the method to which 
this paper is devoted offers a valuable means of training students to appreci- 
ate the meaning of the integral calculus, and moreover a means which is alike 
valuable to students of applied and to students of pure mathematics. It 
quickens interest and develops power. 

The method appears to be due to Duhamel* and it will be referred to in 
this paper as Duhamcl's Theorem. It is contained in the principle {cf. §2 be- 
low) that in the limit of a sum of infinitesimals, any infinitesimal may be re- 
placed by one that differs from it by one of higher order. It was introduced 
into calculus teaching in this country largely through Professor Byerly's, and 
Professor B. O. Peirce's text-booksf and teaching. 

Objections have been raised by Mansion to the ordinary formulation of 
the theorem. J His criticism is chiefly destructive, for while he points out how 
the theorem may be interpreted so as to be correct, the formulation to which he 
is led is, as he himself observes, of no use in practice, and his conclusion is 
that for this reason the only value which the theorem retains is for heuristic 
purposes, — ^as a "principe d'invention." In fact, what Mansion has done in 
this paper is similar to what Seidl§ did in pointing out that a convergent 
series of continuous functions does not necessarily define a continuous func- 
tion, II and in finding out at what point Cauchy's proof of the contrary theorem 
breaks doAvn. But just as Seidl failed to construct from his observations on 
series the principle of uniform convergence, and to introduce that principle as 
a means of investigation, — it was WeirstrassIF who some seven years earlier, in 
two papers not published at the time, had taken this step ; later, through his 
university lectures, the method became generally knoAvn, — so INIansion fails 
to extract from Duhamel's principle as originally stated a theorem fonnulated 

* Mansion (c/.below) refers it to Duhamel, who gives it in his Elements de calcul infini- 
tesimal, Paris, 1856, p. 35, and applies it to numerous problems in geometry and mechanics. 

fByerly, i9!;(reren(i(iJ Calculus, Boston, 1882, ch.x. Infinitesimals; Problems in Differential 
Calculus, Boston, 1895, ch. x B. O. Peirce, Newtonian Potential Function, Isted., Boston, 1886. 

X Mathesis, vol. 8 (1888), p. 149. 

^Abh. der kgl. bayerischen Akad. der Wiss. (math.-phys. CI.), vol. 6 (1848), p. 381. 

jl This fact had already been noted by Abel in his memoir on the binomial series, Journ. 
fur Math., vol. 1 (1826), p. 31C, in the case of the Fourier's series S( — l)''-isin nx/n; and it 
had also been observed by Stokes almost simultaneously with Seidl. 

^ Cf. two papers of the years 1841 and 1842, Werke, vol. 1, p. 67 and p. 75. 
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so simply that it can oasilj he applied to the cases that arise in practice, and 
thus to rescue Duhamel's principle for tlic important place in the calculus 
which it was designed to occupy. To put Duhamel's principle on a firm foun- 
dation and to make clear its importance in the applications of the calculus is 
the object of the present paper. 

2. Duhamel's Theorem. Duhamel gives the following 
Theorem. Let aj + a., + • ■ • + a„ be a sum of positive infinitesimals* 
which approaches a limit when u = oo. Let /Sj + /3., + • • • + y8„ he a second 
sum of positive infinitesimals ivhich differ respectively from the infinitesimals 
of the first sum hi/ infinitesimals of higher order; i. e. let 

lim^ = l. 



n = 3o •*» 



Then the second sum approaches a limit when n = -jo, and this limit is the same 
as that of the first sum : 

lim (^1 + /3. + . . . + /3,.) = lim (ai + a.^ + ■ ■ ■ + a„) . 
To prove the theorem let 



a: 



= 1 + ei, fii = ai + e.ai, 



and let the numerical value of the largest e, be denoted by e : 

I €■ I ^ e «■ — 1 2 . 



Then 



Adding, we get 



Ui — eui ^ /3i ^ oj + eui 
a.2 — ea.2 S yS.j ^ a^ + ea^ 



(l-€)Sa, S 2A ^ (1 +€)Sa,-. 

i = l 1 = 1 1 = 1 



Since lim e = when n = co, each extreme of this double inequality approaches 
as its limit limSa„, and hence 2^„ approaches this limit also, «/. e. d. 



* By an infinitesimal is meant a variable whose limU is zero. 
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3. Applications, (a) Attraction of a Rod. Let it be required to 
find the attraction of a straight rod of variable densitj' on a particle situated in 
the line of the rod. 

Here we start with the physical law of gravitation, namely that two parti- 
cles of masses m and m' attract each other with a force f proportional to their 
masses and inversely proportional to the square of the distance r between 
them : 

. mm' 

Now the rod is not made up of particles, but is a continuous distribution of 
matter of density p. How can we apply the law of gravitation to it? As fol- 
lows : Divide it up into n parts, and denote the length of the i'-th segment 
Xf — Xj_i by AXj. (^Let p'l, p'- denote respectively the minimum and the maxi- 
mum densities in this segment. 7'Then the mass ^Mi of so much of the rod 
obviously lies between p'l^Xi and' />'/ Ax^ : 

p'AXi S A J/, ^ p'JAx,., 

the equality signs holding in the special case of constant density. 

Denote the attraction of the rod which we seek to compute by A, the part 
of the attraction due to the i-th segment by A-4i, so that 

A = 'L AA-. 
t=i 

Now we can approximate to ^Ai as follows : If all the mass AJ/,: were concen- 
trated in a point at the nearer end of the i-th segment, we should obtain an 
attraction KmAMi/xl_i greater than the true attraction A^^, and if we were to 
replace AMf by the larger mass pl'AaCj, the attraction thus obtained would be 
still greater ; hence 

AA, < .^^'- . 

On the other hand, if we were to concentrate the smaller mass /0/Aa;< in a 
point at the further end of the t-th segment, we should obtain an attraction 
Kmp'iAXi/Xi smaller than AJLj. Hence 

Here, m is at the point O: x ~ 0. 
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The expressions that form the extremes of this double inequality have the 
further property that they both differ from AAi by a small percent of AAi 
when AXf is small ; in other AVords, either one differs from AAi by ^.n infini- 
tesimal of a higher order ; for, divide through, say, by the left hand expres- 
sion : 



1 



; ^ j^.i — 

mp'iAXi Pi^v-i 



When ra = 00 and all the Ax,'s approach as their limit, the right hand 
member of this double inequality approaches 1, hence the middle ratio does, 
too, and wq see that Kmp'iAXi/xl differs from A^,- by an infinitesimal of a 
higher order. 

We are now in a position to apply Duhamel's Theorem and thus to evaluate 
the attraction A by integration. Since 

^ = i AAi, 

we may allow n to become infinite, each AAi approaching as its limit, and 
we shall still have 

n 

A = lira 2 A^f. 

n^co isszl 

But here, in dealing with the limit of the sum lAAi, we may replace A^,- 
by any more convenient infinitesimal which differs from it by an infinitesimal 
of higher order. If, now, we can find such an infinitesimal, of the form 
/(Xi)Ax„ then the new limit: 

n 

lim 2 f(Xi)AXi, 

n = 30 i^l 

can be evaluated by integration, for the fundamental theorem of §1 will apply 
to it. 

Such an infinitesimal is suggested by the form of the approximations in 
(1) ; it is : 

mpiAXi 
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where pi is formed for the point x = x,- ; the fuiiction/(x) being 

For, dividing (^1) through by it : 






and the limit of each extreme is 1 . 

Hence by Duhamel's Theorem and the theorem in §1 we infer that 

^ = lim I AA, = lim i .'!^ = r.^^dx, 

and the attraction of the rod is formulated as a definite integral. If, in particular, 

/) is a constant, then 

fdx (\ 1\ 

The mass J/ of the rod is here p(b — a), and thus 

m3f 



A = K 



ah 



The rod attracts, therefore, like a particle of the same mass concentrated, not 
at the centre of gravity of the rod, but at a distance from O which is a mean pro- 
portional detween a and b. 

(b) Attraction of Any Material Body. Let us now determine the attrac- 
tion of an arbitrary material body on a particle. We begin by analyzing the 
method. 

Let a system of Cartesian coordinates be assumed, and let the coordinates 
of the particle m be denoted by (a, b, c) . We wish to find the components 
X, Y, Z along these axes of the attraction which the body exerts on m. Begin 
with the component X. Divide the volume V occupied by the body into n 
small pieces A F^- . . . A F„ and denote the component attraction of the mass 
AJ/f contained in the volume A V,- by AX^. Then 

n 

A'= 2 AA-;. 

(■ = 1 
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The method is now this. If we allow n to become infinite and divide F"up in 
any manner such that the maximum diameter of each A Vi approaches as 
its limit, then, of course, 

A"= iim i A A';-. (2) 

n = oo i = l 

This limit, as it stands, we do not know how to compute ; but we do know how 
to compute the limit of a sura of the following form : 

limif(Xi,yi,Zi)AVi, (3) 

wherey" (^cc, y, z) denotes a function continuous throughout V. For, by the 
fundamental theorem of the integral calculus* referred to in §1, this limit exists 
and is the triple integral off{x,y, z) extended throughout V: 



III 



f(x,y,z)dV. 

V 

It may be computed by means of three successive simple integrations ; for 
example, by the formvila, 



rb ry^ rzt 

dx dy I f(x, y, z)dz. 

Ja Jiio Jz, 



If, then, we can succeed in throwing the limit (2) by means of Duhamel's 
Theorem into the form of the limit (3), then we can evaluate (2) by simple 
integration. This can, in fact, be done, and we turn now to the proof that it 
can be. Assume for the present that the body lies wholly to the right of the 
plane a; = ; i. e. that at no point of the body x is negative. 

We shut AAi in between two expressions differing each from AA^- by a small 
percent of AA",. The form of these approximations will suggest the form of the 
function f{x, y, z) to be employed. Let p'j, rj', a'j denote respectively the 
maximum values of the density p in AF^, the distance from (a, b, c) to a point 
of A Fi, and the angle a which a ray ditiwn from (a, b, c) to a point of A Vi makes 
with the positive axis of x ; and let p'l, ?•;, o; denote respectively the minimum 
values of these functions. Then it is clear that if a particle of mass Ailifj were 
concentrated at the shortest distance r'l from m to a point of A Vi, and were 

* Of., for example, Goursat, ch. 7. 
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placed, not necessarily at a point of AF'j, but on ai-ay making the minimum 
angle a'f with the x-axis, the component attraction 

K ;;— cosa;- 

' t 

:vould be greater than AJ^^ ; and if AMi were replaced by the larger value 
Pi' A Vi — we have as before 

p'iAVi^AMi^p'JAVi — 

then the inequality would only be strengthened. In a similar manner we 
obtain an expression that is less than AX^, and hence the double inequality 

mp'iAVi ,, . _^ mp'i'AVi , ,,, 

K ^' , ' cosaj' < AXi < K-f-^, — ? cosa^ (4) 

ft **» 

The first and last members above suggest a form of the summand in (3) 
such as is desired, namely, the summand obtained by setting 

f(x, y, 2) = /c-^cosa. 

For, dividing (4) through by this function formed for any arbitrary point of 
A Vi and multiplied by A F^, we get 

p\ i\ cos a'l AXi Pi »1 cos aj 



< 



PiT'i cosai ' ^rnp^^^^^_ p^rfcosa* 



< "':;• • (5) 



n 



When A Vi shrinks down toward a point as its limit, each extreme of this double 

inequality approaches 1 as its limit, and hence the middle term must do so also. 

Thus the conditions of Duhamel's Theorem, §2, are satisfied and we have 

J y (6) 

f f f P(^ — a)dxdijdz 



Km 



V 

with similar fornmlas for I^and Z. 

If the solid lies wholly to the left of the plane x = 0, the component at- 
traction in the direction of the negative axis of x can be computed in precisely 
the same manner ; and if the solid lies on both sides of the plane, then the 
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attraction can be computed for each piece separately and the results, affected 
by the proper algebraic signs, added. The final formula will be the same as 
the above, (6). 

4. Critique of Duharael's Theorera and its Application in §3. 

The infinitesimals a^, /3,- that appear in Duhamel's Theorem depend on two 
variables, i and n, and the hypothesis 

V ^' 1 
lim — = 1 



H=ao *"» 



tacitly assumes that i varies with n according to some specilied law, or per- 
haps according to any law whatever ; it is not explicit on this point. The proof 
tacitly assumes that lim e = when n = oo . 

An example will make clear questions which a critical examination of the 
law and its proof suggest. * Consider the curve 

y = 1 + n*xe-"'*' (7) 

in the interval S .x ^ 1. The graph of the curve is as indicated in the ac- 
companying figure. For any specified value of x in the interval the ordinate 
of the curve approaches 1 as its limit when n = ao ; but the 
peak rises higher with increasing values of n and moves at the 
same time further and further to the left toward the axis of 
ordinates. At first sight it might appear as if the area under 
the curve necessarily approached as its limit the area under 
the limiting curve, y = 1; i. e. the area 1. But another factor 
comes in, namely, the ever increasing height of the peak, and 
so we see that it is not certain that the limit of the area under the curve is 
what it seemed to be. In fact, let us determine it ; it is : 




A„= i ydx = (1 + n^xe-'''''')dx 
Jo Jo 

= (x - i e-»'-)|; = 1 -H i(l - e-"'), 



and lim A^ is seen to be IJ instead of 1. 

Let us now divide the interval (0, 1) into « equal parts and through the 
points of division draw lines pai-allel to the axis of ordinates, thus dividing the 

* This is the geometric form of an example given by Mansion, I. c. 
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area under the curve (7) into strips. As n increases, some of these strips 
will run up into the peak of the curve, so that their maximum altitude will in- 
crease indefinitely with n. In fact, let A.t,- = 1/n, Xi = i Ax^. Then 

yj = 1 + nie-'' 

and if i is held fast while n increases, y,- becomes infinite with n. 

As the infinitesimal a^ of Dahamel's Theorem we will take the area of the 
t-th one of the above strips, namel}'. 



a,- = / (1 + n^xe-»'''')dx 



as /Sf we will take the part of the area of this strip lying below the limiting 
curve y = 1, namely, 

A- = Ax,. 

Then if we fix our attention on the strip that contains a given point, as for 
example, the point x = X, y = 0, and allow n to become infinite, i will become 
infinite at the same time and x,, Xi_-^ will both approach X as their limit. 
Hence we shall have 

Urn ^* = lim - "" 



n=^ ai „=« 1 ^ ^(e-'''*J-. - e-»''A 



and the conditions of Duhamel's Theorem, as they might readily be understood, 
are fulfilled. But the theorem would then be false, for 

n n n 

S «, = S Ax, + S i[(l - e-i) + (e-i - e-*)+ • ■ • +(e-(»-i)^ - e""')] 
»=i «'=i »=i 



and lim 2 o,- = 1^. 

n s= 00 i = I 

The explanation is simple. For the case in hand the maximum e^, taken 
numerically, namely e, does not approach zero. Or again, it is possible to al- 
low i so to vary with n that /i^,/a,- does not approach 1 ; for if we hold i fast, and 

express Ax,- as - , then — = [1 + :^ (e-(»~i5' — e-*')]-^ and ^i/aj approaches 
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as its limit when n becomes infinite. Here, the number of the strip 
remains constant and the strip moves further and further to the left.* 

It is clear from this example that the formulation and the proof of the 
theorem leave someth-ing to be desired. t 

The foregoing critique applies to Duhamel's Theorem itself. Let us noM' 
turn to its application in §3. First, the restriction that all the infinitesimals 
must be of like sign made a clumsy division into special cases necessary, the 
general result being identical with that obtained in each of the special cases. 

A more serious difficulty is the following. Suppose the body extends 
up to the plane x — 0. Then for a AFj abutting on this plane a(' = j^7r and 
cos a'- = 0. Hence it is not true in this case that the limit of the left hand 
member of (5) is unity, and the hypotheses of Duhamel's Theorem are not ful- 
filled. It would be possible, to be sure, to take the boundary of the solid 
distinct from the plane a; = 0, but near to it, and to allow the boundary sub- 
sequently to approach this plane as a limit- But to be obliged to take such 
precautions as these and go through with the consideration of a passage to 
a limit would detract very largely from the simplicity which is desirable in 
the case of so important a theorem. For the reasons here given we modify the 
statement of Duhamel's theorem in such a manner that it will apply to the 
cases that arise in practice and at the same time retain a satisfactory degree 
of simplicity. 

5. Revised Formulation of Duhamel's Theorem. The impor- 
tant applications of Duhamel's Theorem in practice are to cases in which we 
wish to show that a given limit of a sum : 

n 

lim 2 o,- 

n=:oo f=l 

• This example is closely related to the question of uniform convergence and the integration 
of a series term by term ; ef. the author's geometric treatment of this subject, Bulletin Amer. 
Math. Soc, ser. 2, vol. 3 (1896), p. 59. 

fThe formulation may be made precise and the theorem correct by requiring that 
lim ft/a, = 1, no matter how i varies loUh n, as n becomes infinite. For it can then be readily shown 
that lim «=0, and thus the proof of §2 will hold. But this formulation is of little use for prac- 

tical purposes, since there is no simple sufficient condition which insures that the hypotheses 
of the theorem are fulfilled ; and even if such a test were found, the theorem would still be in- 
applicable to simple cases that arise in practice, as will presently appear. The formulation 
given below (§5), on the other hand, is not open to this objection. 
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is equal to a definite integral,* either to 

n rb 

lim 2 f{Xi)Axi = / f{x) dx, 

n = 00 f = I Jo. 



or to 

lim 2 f{x„y^)AS, = [ ff{x,y)dS, 

n=: CD i=:l J J 

s 

or to 

lim if{Xi,i/„z,)AVi= f f ff{x,y,z)dV. 

V 

We will restrict ourselves for simplicity to the first case. 
Definition, Let fj be defined by the equation 

a,: =/(a:.-)AXi + ^lAXf. 

Then the infinitesimal o^ is said to differ uniformly from the infinitesimal 
f(Xi) AXi by an infinitesimal of higher order than Axi if, a positive quantity e 
being chosen arbitrarily small, it is then possible to find a second positive con- 
stant B such that, no matter how the interval a ^ x ^ 6 be divided into n parts, 
we always have 

provided merely that the division has proceeded far enough so that each 
AXi < B. Here a^ may be positive, negative, or zero.f In the case of a double 
or triple integral ?,- will be defined by the equation 

«i =/(a;» yd^^i + tASi or a,- =f(Xi, y,-, Zi)A Ff + fiA Fj, 

where ASi, A Fi denote respectively a piece of the surface or the volmne into 
which the i-egion of integration has been divided. 

*The notation here is not the same as in §2, the present o,- corresponding to the former Pi 
and the present /(x,)Aa^ to the former a,-. 

flf one wishes to admit negative Axj's to the consideration, then the last relation must read 

I AXi I < «. 
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Theorem. Let 
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«! + «2 + • • • + a„ (^) 

be a sum of infinitesimals and let a,- differ uniformly by an infinitesimal of 
higher order than ^Xifrom the summandf(Xi)AXi of the definite integral 



i: 



f{x)dz 



{B) 



of the function f {pi), this function being continuous throughout the interval 
a ^x^b. Then the sum (A) approaches a limit when n = co, and the value 
of this limit is the definite integral (B) : 



Prooj. Since 
and 
it follows that 



» rb 

lim 2 Oi = I f{x)dx. 

\Ki\ < e, AXf < S, 



n » n 

2 a; = 2 f{Xi)AXi + 2 ^iAxi, 

t = l 1=1 i = l 



|2ai - I.f(Xi)AXi\ S 2|?i|Ax,- < 2 eAx.- = €(6 - a), 



isrl 



i=l 



1 = 1 



i = l 



(8) 



when AXi < 8. Now the sum 1f(Xi)AXi converges toward the integral (B) 
when^ = 00, and so 



2 f{Xi)AXi - / f{x)dx 



»=i 



< e. 



(9) 



when AXi < B'. If 8 is chosen at least as small as 8', relations (8) and (9) will 
hold simultaneously as soon as AXi < B, t = 1, • • • n. From (8) and (9) 
we infer that 



n rb 

• = 1 Ja 



x)dx 



< €{b — a) + e, 



when AXi < B, and this relation is coextensive with the assertion of the theorem 
we wished to prove. 

Both theorem and proof apply without modification to multiple integrals. 



174 OSGOOD [July 

6. Application of Duhamel's Revised Theorem to the Ex- 
ample of §3. The formulation of the problem down to the point of the 
double inequality (4) is the same as in §3. We now identify the infinitesi- 
mal A^- with the a^ of §5 and set 

f{x,y,z) = K-^cona. 

The functions p, r, cos a are continuous throughout the Avhole volume F in- 
clusive of the boundary, and r never vanishes in this region. Hence p and 
r-^ cos a are both uniformly continuous throughout V, and the oscillation of 
each in A Vi can be made uniformly small by choosing all the A V/s sufficient- 
ly small. If, then, 17 denote an arbitrarily small positive quantity, the oscil- 
lation in each A Vi will be less numerically than 7; for all divisions of V for 
which the maximum diameter of any AF^ is less than a suitably chosen posi- 
tive constant S. Moreover the maximum numerical value G of either function 
in V is finite : 

\P\^G, 



Thus we shall have 



^ G. 



p'l = Pi + o?\ 

COSaS _ cosa; 
r? - A ^ " 
where pi, r,-, coso^ apply to an arbitrarily chosen point of AF,- and 

\e?\<v, \er\<v- 



Hence 



mp'/A Vi , mpjA Vi , ^„ . ,^ 

<^ \.,, cosa^ = K ^'^ cos% + C?A1',-, 



where 

^V = Km{ri-^ cosai^(" + /j.-^f + 62^6^^) and | r/| < Km{2G + 7;)^ = e. 

A similar relation holds for the left hand member of (4) . Hence if we 
subtract the expression 

m/3jA Vi 
K—!—„ cosaf 

from each member of (4), the new extremes will each be less numerically than 
eA Vi, and we get 

I AAi — K —!-^ — • cos«i I < eA Vj. 
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This relation expresses precisely the fact that 

K 5 COSai 

differs uniformly from AJ^j by an infinitesimal of higher order. 

7. The Axioms of Physics. Gravitation. We have chosen for pur- 
poses of illustration a problem in gravitation. The method set forth applies 
in the same manner to problems in centres of gravity, moments of inertia, fluid 
pressures, kinetic energy of a rigid body (or, generally, of any continuous dis- 
tribution of matter), work done on a rigid body by gravitational forces, etc. 
It is interesting to inquire in each case what the physical laws assumed are, — 
what is definition and what is proof. 

Consider the problem of §3. Here, to begin with, the law of gravitation 
for two particles is assumed. But is that all? Have we, with the aid of this 
law alone, computed the attraction of the given body? By no means. We 
have assumed as further axioms 

(a) that a resultant attraction for the whole body and for any piece of 
it, exists ; 

(6) that when the body is divided up in any way the resultant attraction 
due to the attractions of all the component pieces is equal to the resultant at- 
traction of the whole body ; 

(c) that the intensity of the resultant attraction for a given piece is less 
than the force which a particle of the same mass as this piece would exert if 
placed at the minimum distance from the given particle to a point of the piece, 
and greater than if that particle were placed at the maximum distance ; 

(d) that the direction of the resultant attraction of any piece lies within 
any convex cone that can be drawn with the given particle as vertex, contain- 
ing the piece in question.* 

Moment of Inertia. An analogous set of axioms is employed in finding 
centres of gravity of continuous distributions of matter, and moments of inertia 
may be treated in the same way. Another point of view regarding moments 
of inertia is as follows. Consider, for example, the moment of inertia of a 
plane lamina of variable density about an axis perpendicular to its plane and 



♦This axiom regarding the direction of the resulting force is not the same in form as the 
one employed In §3 ; but it is evident that the inequality (4) of that paragraph can be deduced 
from it. 
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piercing the plane in the point 0. Here we begin with the definition of the 
moment of inertia of a particle of mass 7n situated in the plane, about O as 

where r denotes the distance of the particle from O, and we define the moment 
of inertia of a system of such particles as the sum of the moments of inertia of 
the individual particles. Now instead of considering the moment of inertia 
of the lamina about as a fundamental concept and proceeding to compute its 
value as above suggested, we may regard this moment of inertia as a quantity 
to be defined and for the purpose of the definition lay down requirements sug- 
gested by the physical meaning of the moment of inertia of a particle or system 
of particles. These requirements we may state as follows: The moment of 
inertia of the lamina, and of any piece of it, shall be so defined, 

(a) that when the lamina is divided up in any way, the sum of the mo- 
ments of inertia of the component pieces shall equal the moment of inertia of 
the whole lamina ; 

(6) that the moment of inertia of a piece shall lie between the moments 
of inertia of a particle of the same mass when concentrated at the nearest point 
of the piece and at the farthest point respectively. 

The question that now presents itself is whether such a definition is pos- 
sible, and if so, how the moment of inertia may be computed. Let us see. 
Suppose that the definition is possible. Denote the moment of inertia of the 
whole lamina hy I, that of one of the component parts into which it is divided 
by Ai^-. Then, if Aif/j-, A/S*;, p'l, p", r[, r'l denote respectively the mass, the 
area, the minimum and the maximum densities, and the minimum and the 
maximum distances from to a point of A/iS',-, we shall have 

/ = I A/i = lim 2 A/j, 

pjA^rf < Aii < f^!t.S,r^t. 

The extremes of the last inequality suggest the function 

/(x, y) = pr'' 

formed for an arbitrary point (x,-, y^) of the region AxSi. Now p and r are both 
continuous functions of x, y, throughout the whole region of the lamina ; hence 
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they are uniformly continuous, and thus \ve are able to show, just as in §6, 
that A/j differs uniformly from 

by an infinitesimal of higher order. All of the conditions of the revised form 
of Duhamel's Theorem, §5, are fulfilled regarding the infinitesimals Aj/and we 
see that /must have the value 

/=liiu 2 pir'iASi= f fpr'^dS. 

s 
This is a necessary condition for the existence of a moment of inertia which 
shall satisfy the given requirements. Conversely, it is sufficient, as is readily 
shown by the law of the mean. Hence we may define the moment of inertia 
of the lamina as 



■=IJp^<i^. 



S 
and this shall be its definition. 

Length of Arc of a Curve. For the elementary student the length of 
the arc of a curve is a fundamental concept. He is ready to accept the axiom 
(i.e. physical law) that the length of an arc of a plane convex curve, when 
the arc lies wholly on one side of the straight line which passes through its 
extremities, is less than the length of a broken line enveloping it and having 
the same extremities ; while on the other hand, the length of the arc is greater 
than that of the chord . The ordinary formula for the length of a n arc of a 
plane curve may then be proved as in §6, On the other hand, for the pro- 
fessional mathematician the length of the arc of a curve is a matter of defini- 
tion, and the requirements on which the definition depends lead to considera- 
tions analogous to those above discussed in the case of the definition of the 
moment of inertia. 

8. Conclusion. In closing, one word about the place of these theorems 
in elementary instruction in the calculus. It is advisable that the integral as 
the limit of a sum be introduced at an early stage and illustrated by numerous 
problems in the determination of volumes of geometric solids. The solid in 
question is sliced up and the actual slice is compared with an auxiliary slice 
whose volume is readily computed in the form 

f(Xi)AXi 
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and which is so chosen that the total volume of these auxiliary slices, 

I /(cc^AXi, (10) 

obviousl}' converges toward the volume sought as a limit. But the limit of 
(10) can be determined by integration, §1. In many cases, as for example in 
finding the volume of a regular pyramid or of an anchor ring, a formal proof 
may readily be given, and it is desirable that some such proofs be given. But 
the observational faculties of the student should also be exercised and formal 
proofs, when clumsy, may well be waived in favor of greater stress on the in- 
tuitional side of the problem. 

Proceeding now to the next stage, the treatment of physical problems, we 
find that intuition is less clear. The student may be willing to grant the for- 
mulation of the attraction of the rod in § 2 as the definite integral there found 
for it, but his reasoning is more likely to be based on formal analogy, such as 
led him a year earlier to think that sin 2x is equal to 2 sin x, than on an estimate 
of the relative magnitudes of the infinitesimals with which he is dealing. It is 
here that Duhamel's Theorem in its original form (§2; gives him a rigorous 
principle by means of which to determine whether his processes are correct. 
I say a rigorous principle, for the principle meets the highest demands 
for rigor of which the student is capable at this stage. What is the state of 
affairs at present regarding the treatment of these physical problems ? The in- 
finitesimal is treated as an inconceivably small constant, and an infinitesimal 
which does not succeed in making itself persona grata is neglected. Shall 
we hesitate to replace this barbarous method by one right in principle and 
easy of application simply because the whole truth can not be told at this 
point ? The systematic treatment of physical problems by the method of limits 
and Duhamel's Principle means a step in advance in calculus teaching, and it 
is possible and desirable for this generation of mathematicians to take this 
step, both in college courses and in the courses of our technical schools. 

Finally the student of mathematical physics and of higher analysis, who 
has reached the stage at which questions of uniform convergence are treated, 
will find in the revised form of Duhamel's Theorem a method rigorous accord- 
ing to the highest standards of today. 

Harvakd Untvbrsity, 

Cambmdge, Massachusetts. 



